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We study the phase diagram of a chiral random matrix model with three quark flavors at finite 
temperature and chemical potential, taking the chiral and diquark condensates as independent order 
parameters. Fixing the ratio of the coupling strengths in the quark-antiquark and quark-quark 
channels applying the Fierz transformation, we find that the color-flavor locked (CFL) phase is 
realized at large chemical potential, while the ordinary chirally-broken phase appears in the region 
with small chemical potential. We investigate responses of the phases by changing small quark 
masses in the cases with three equal-mass flavors and with 2+1 flavors. In the case with three 
equal-mass flavors, we find that the finite masses make the CFL phase transition line move to the 
higher density region. In the case with 2+1 flavors, we find the two-flavor color superconducting 
phase at the medium density region as a result of the finite asymmetry between the flavors, as well 
as the CFL phase at higher density region. 
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I. INTRODUCTION 

Mapping of QCD phase diagram at finite temperature and densityfTj-Q is one of the most challenging issues in the 
theoretical and experimental physics and is significant to the heavy-ion collision experiments and the structures of 
the neutron stars. 

At finite temperature and zero or small baryon density, a number of investigations on the QCD phase transition are 
made both with the lattice QCD simulations @ and with the model calculations (if. Lattice QCD simulations suggest 
that the phase transition becomes smooth crossover in the realistic case with two light and one heavier quark flavors, 
and many model calculations are consistent with this result. At finite density, however, the situation is uncertain 
because lattice QCD simulations are still challenging at finite chemical potentials with low temperatures Q- In such a 
region, it is important to employ models for qualitative and quantitative calculations on the phase diagram. Naively, 
a large baryon number density may cause overlaps of baryons, which invalidates a concept of confined colors in a 
baryon, resulting in dcconfinemcnt, which may be followed by chiral phase transition. 

Furthermore, at asymptotically large density, the appearance of a color-superconducting (CSC) phase is also ex- 
pected, where a weak coupling theory is applied and the Cooper instabiltiy of the Ferimi sphere is inevitable 0, 
The color superconducting phases are characterized by the diquark condensates. One-gluon exchange interactions tell 
us that the color- antitriplct channel is attractive. Therefore, with the Pauli principle, condensates in the color- and 
flavor-antitriplet and spin antisymmetric channel is expected . Such condensates can be expressed as 

SAA> = (l^ c C^ 5 T A X A >lfj) , (1) 

where ip is the quark field, and ip c = i[> T C with C the charge conjugation operator, ta and A,4, where A = 2, 5 and 
7, are the antisymmetric generators of SV(Nf) flavor and SU(A r c ) color groups, respectively. 

One of the most striking features in the CSC phases is the formation of the color-flavor locking (CFL) condensate. 
At sufficiently high density, the finite current quark masses for up, down and strange quark flavor are neglected so 
that the system can be treated as the chiral limit. In the CFL phase, which is characterized by saa' = A5aa' with 
nonzero A, SU(3)l® SU(3),r(8> SU(3) c symmetry of the system breaks down to its subgroup of SU(3)l+_r+c- This 
breaking pattern is possible owing to the miraculous matching of the (effective) number of flavors JVy = 3 and the 
number of colors N c = 3. 

To investigate the QCD phase diagram at finite temperature and chemical potential, chiral random matrix (ChRM) 
models provide us of a qualitative way from a viewpoint of the symmetry [lol [Tlj . In a conventional ChRM model, 
the Dirac operator is set as a random matrix, which has the same symmetry with QCD, and the partition function 
is defined as the average of the determinant of the Dirac operator over the matrix elements with the Gaussian 
distribution. The random distribution of the matrix elements mimics the complex dynamics of the gluon fields. The 
Gaussian model can be solved exactly in the thermodynamic limit. Although the model is constructed in such a simple 
way, the resulting phase diagram has a rich structure. In the chiral limit, the phase transition becomes second-order 
in the small chemical potential \x region, while it becomes first-order in the large fi region. First- and second-order 
phase transition lines are connected at the tricritical point (TCP). This result is consistent with NJL models (l2l. fl3|. 

The extension of the ChRM model to the case with the CSC phase has already been studied by Vandcrhcyden 
and Jackson 140. In their study, the Dirac matrix is extended to have the indices of color and spin explicitly with 



the real random matrices corresponding to real gluon fields in QCD. After the integration over the random matrices, 
the model produces the quark-quark interaction terms, which are responsible for the diquark condensates, as well as 
quark-antiquark interaction terms responsible for the chiral condensates. The resulting phase diagram has the diquark 
condensed phase at large /i region, while the chirally-broken phase at small /x region, if the ratio of the quark-antiquark 
and the quark-quark coupling is taken so that the Dirac operator of the model has the same symmetry with that 
of QCD. Note that because the model in Ref. [l4[ contains two quark flavors, the CSC phase is the two-flavor color 
superconducting (2SC) phase, where only two of three colors of fermions participate in the diquark pairing. 

It is then natural to ask whether it is possible to extend the ChRM model to the case with three flavors, and 
whether the CFL phase can appear as the ground state in a high density region. We answer "yes" to this question by 
constructing the ChRM model containing three flavors and colors, and show a phase diagram with the chirally-broken 
and the CFL phase. As a simple application of this model, we also focus on the response of the model by changing 
the quark masses. By setting the strange quark mass different from the other two quark flavors, we observe that the 
2SC phase appears on the phase diagram at the moderate values of the chemical potential, as well as the CFL phase 
in the larger chemical potential region. 



1 Also, there are studies on the diquark condensates with N c = 2 using ChRM models [Hill. 
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This paper is organized as follows. We introduce an extended ChRM model with chiral and CFL condensations 
in sec. 2, and derive its effective potential in sec. 3. The model phase diagrams are presented and discussed in the 
case with three equal-mass flavors and 2+1 Savors in sec. 4 and 5, respectively. Sec. 6 is devoted to a summary and 
discussion. 



II. RANDOM MATRIX MODEL WITH CHIRAL AND DIQUARK CONDENSATIONS 

In this section, we introduce a chiral random matrix model which mimics QCD partition function with three quark 
flavors, extending the two- flavor case in Ref. [l^. We denote three quark masses by m/ with / =u, d and s. 

Keeping in mind the (extended) Banks-Casher relations, which relate Dirac soft modes not only with the chiral 



condensates [1 5} . but also with the diquark condensates (16j. we consider the truncated Dirac matrix D in low- lying 
quark excitation, or zero-mode space. We assume D can be separated as D = R + C, where R is a random part, 
which represents the complex gluon dynamics, and C the non-random, deterministic part, which is responsible for 
the matter effects. 

For simplicity, we first set the matter effects turned off, i.e. C = 0, and focus on the random matrix R. In this case, 
the truncated Dirac matrix should have the symmetries of the Dirac operator in the vacuum, the chiral symmetry, 
{75, R} = 0, and the anti-hermiticity, R' = —R. Within these restrictions, the Dirac matrix generally has nonzero 
matrix elements only in the off-diagonal blocks 

in the chiral representation, 75 = diag(+l, —1), where W is a complex matri^fl. In the conventional ChRM mo dels p^j. 
W is taken to be a general complex matrix whose elements are independently distributed according to the Gaussian 
distribution. 

To consider the diquark condensations, however, it is crucial to treat the color and spin indices explicitly. Following 
the construction in Ref. [l4| . we express W as a direct product of the spin, color and zero- mode matrices, whose total 
dimension is 2 x N c x N, where 2 is the size of the spin space, N c the color space, and iV the zero-mode space. We 
adopt the form of W as 

W =A va (a„®\ a ), (3) 

where a v = (1,—iai) with Ui the Pauli matrix, A a is a generator of SU(iV c ), and A va is a N x N random matrix. 
Since the random matrix A va corresponds to the gauge field in QCD, we choose A ua to be a real matrix. 

The matter effects are introduced as the non-random external fields in the Dirac operator. A simple way (Til] is to 
add a constant matrix C to the random matrix @ with 

C=(£X)=«®Tb, (4) 
where the 2 x iV r x iV-dimensional matrix lu is defined as 



(H + iT)l N/2 ® lspin® 1jv c 

(/J,- iT)l N/2 <8> lspin ® 1jv< 



(5) 



T and \i are an effective temperature and quark chemical potential, respectively. l sp i n denotes 2x2 identity matrix 
in the spin space. A total Dirac matrix, D = R + C, also has chiral symmetry, {£), 75} = 0, but not anti-hermiticity, 
D' 7^ — D if /1 7^ 0. Two relative signs between T and /i corresponds to the two lowest Matsubara frequency, ±7rT. The 
inclusion of two signs reproduces the invariance of the partition function under the charge conjugation transformation, 
(i -> -fi. 

Using the Dirac matrix D, the ChRM model partition function is defined as 

Z= [dA] [] det(D + m f )e- 2Ns2 ^,J A "?)\ (6) 
J f=i 



2 One can take W generally to be an JV+ X N— rectangular matrix. In this case, the Dirac matrix has \v\ = |JV+ — iV_| exact zero 
eigenvalues, which represents the index theorem with the background gauge field having the topological charge v. Exploiting this fact, 
we can introduce the effect of the axial anomaly in the ChRM models [l^, [2(1 . In this study, however, we always take W to be square 
and the anomaly effect is neglected. Sec the discussion in sec. 6. 
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where the integral is defined over real elements of random matrices A va with the Gaussian weight. The parameter E, 
which fixes variance of the Gaussian distribution, gives a scale to the model. Generally, X may change for each a and 
is, but by ensuring the color and Loerntz symmetry, their values should be equal. 

Before solving the model, we make two remarks on the treatment of the ChRM model comparing to that in Ref. [l4[ ■ 
First, in Ref. 14|, the authors examine not only the form of the random Dirac operator (|2|), but also the case where 
the Dirac operator breaks the symmetry which QCD Lagrangian holds. In such general cases, the random matrix R 
is taken as 



(7) 



Ca 



where the chiral indices a, (3 = R, L, 
dimension, C = 1, . . . 16. A set of T 



vCa 

A af3 



is a random matrix, and Tq is the independent gamma matrix in four 
can be separated into the subsets forming a Lorentz scalar, pseudo-scalar, 



vector, axial- vector, and tensor. In Eq.@ and ([3]), we have chosen X^S to be nonzero only for the vector content of 
the gamma matrices. This is a natural choice because we consider that the random matrices model the gluon fields, 
which form a Lorentz vector. Indeed, nonzero components for the scalar, pseudo-scalar and tensor break the chiral 
symmetry explicitly, and that for the axial- vector does the antihermiticity. If such non QCD-likc random matrices are 
allowed, one can vary the ratio of the coefficients of the quark- ant iquark and the quark-quark interaction channels, 
which is denoted as B/A in the next section. The evolution of the phase diagram with B/A changed was completely 
investigated in Ref. [14| for the case with two flavors. In our study, however, we only focus on the case in which the 
model has the same type of the interaction with QCD, because the condition fixes the topology of the phase diagram 
without ambiguity and we can solely concentrate on the response of the phase diagram by changing the quark masses. 



Second, the scheme of the temperature dependence is different from that in Ref. [14| . We use the matter effect matrix 
proportional to the identity in the flavor space, while, in Ref. [Til ], the sign of the temperature T is antisymmetric in 
the space of two flavors. As a result of this difference, our model can not be reduced to the model in Ref. [l4| in the 
two flavor limit. Further discussion will be given by comparing the effective potentials in the later section. 

Finally, note that the partition function (|6]) has SU (Nc) C <8>SU (N f) z,<g)SU (N f) l global symmetry, but not local 
gauge symmetry. To be exact, it is then appropriate to call the diquark condensed phase a BEC state, not a BCS 
state. In the remaining part of this paper, we use the word, "diquark condensates" to indicate such condensates, but 
discuss them comparing to the BCS states expected in QCD at finite density. 



III. THE EFFECTIVE POTENTIAL 

In this section, we derive the effective potential of the ChRM model defined in Eq.©. The derivation is almost 
parallel to that in Ref.[l4|. We present the derivation in three steps, and then make a few remarks. 

A. Gaussian integral 

The first step is to integrate out the Gaussian integral variables A a ^ . For this purpose, we first express the 
determinant in the partition function ^ in the form with the fermion integrals: 



f r 

J det(Z) + m f )= / [<# + ][#] exp 
dip^ldip] exp 



/=i 



irlAt'-^'iC + rrif)^ 
! 



(8) 



where ipf = (ipR, ip{j T and ipf = (ip{} , ip^) are 4 x N c x N Grassmann vectors, and 



(9) 



is a fermion bilinear. 
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Applying the Gussian integral formula up to the constant, J dxe ax +l3x = cxp(/3 2 /(4a)) , to the A^ v integral 
separately for each indices i,j, a and v, we obtain analytically 



J\dA\e- u ^"e~ 2N ^ ( ^) 2 = exp 



1 



87VE 2 

where the summations over a and v should be understood on the right-hand side. 



(10) 



B. Fierz transformation 



In this step, we expand the square of the fermion bilinear ( J l J v ) 2 and realign the four point vertices, by applying 
the Fierz transformation formula. The square of J l ^ y is expanded as 

+ 4\a v \ a ^ Lj r L \<T»x a r Lj + ^4^i 3 rJA^r Rr (ii) 

The first term represents the quark-antiquark interactions, and the other terms the quark-quark interactions. The 
former is responsible for the formation of the chiral condensates and the latter for the diquark condensates. 

Using Fierz transformations, these four-fermion terms are rearranged so that in each fermion bilinear terms, zero- 
mode indices i and j are contracted. At this point, we assume that the chiral condensates are formed only in the 
color-singlet, scalar channels, and that the diquark condensates are formed only in the spin-antisymmetric, flavor- 
and color-antitriplet, scalar channels. With these assumptions, relevant interaction terms are drastically reduced and 
shown explicitly a^H 

^nW^aKM^^l = - 2{N l7 1] < f ^ 9 M + ■■■ ( 12 ) 

for quark-antiquark channels and 

for quark-quark channels, with the same except L^R, where the charge conjugation fields are defined as tp c = C^ T = 
(?/)£, V , |j)' r an d V' c = ip T C = {i ,C R j V'l')- The dots denote the terms irrelevant to the formation of the condensates we 
focus. 

Neglecting the irrelevant terms, we finally obtain the simple form of four-point interaction as 

- G a ^ l t a X a ^I^t a X a ^ l + {L -> Pi}, (14) 
where we have defined coefficients G x = ^ and Ga = ■ 



C. Bosonization 

We apply the Hubbard-Stratonovich transformation formula, e^ 1 ^ 2 / 4 " = J d ze - a \ z \ 2 +0i z +02 z ; to the rearranged 
four-point interaction (| 14[) . For simplicity, we make further, but moderate assumptions in the formation of the chiral 
and diquark condensates. For chiral condensates, we assume that only the flavor-singlet condensates are formed, and 

for diquark condensates, that only the color-flavor-locked condensates are formed, i.e., s AA ' = \4''l t aX A '4 i l) oc 5 aa >. 



3 For the Fierz transformation formulae, see, for example, Ref. [f| . 
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These assumptions allow us to bosonize the fcrmion vertex ((10]) as 
1 



exp 



8NY 2 



/{ NY 2 NY 2 
[d0][dA]esp (-—2^-^(1. 



2G A 



L \2 
A\ 



lA^I 2 ) 



x exp 



{ff^ f M + l£l>k<t>f) - \ (a^taXatPl + ^ L r A XAr L K) - {L 



R} 



I 



[dcf>] [dA] exp 



NY 2 



2G X 2Ga 



|A*| a 



IA21 2 ' 



) exp 



(15) 



where the measure [d</>] [dA] = !!/=« d s 11.4=2 5 7 ci</!>/<iA^dA^. We have defined the Nambu-Gorkov spinors by 




and 



and the 2 x TVc x Nf(= 18)-dimensional order parameter matrix 5 by 

5 = 



(j>l Nc A\t a \a 



A^taXa 



(16) 



(17) 



(18) 



where <j) = diag(0 u , <pd, 4>s) is a matrix in the flavor space. 

By representing the mass and matter effect terms in the Nambu-Gorkov basis, we obtain the partition function as 



Z 



[# T ][#][#pA] exp 



NY' 
2G~3 



2|<M 



NY 2 
2Ga" 



(|A 



L |2 



| A H,2- 



) exp 



S + M C 
C S^+M^ 



where G = diag(w, —lo) and M is the extended mass matrix 

f fnl Nc -i] A t a X A 
r) A T A \ A ml Nc 



M = 



(19) 



(20) 



with mass matrix rh ~ diag(m u , m^, m s ) in the flavor space. We have introduced the external field r] A , which should 
be zero in the end of the calculation. It is useful to clarify the meaning of the order parameters. 
Finally, the evaluation of the fermion integral is straightforward, which yields 



Z = 



(NY 2 NY 2 
-^2|^| 2 -^-(|Ai| 2 

[d(b}[dA}e- 2NN f N ^ A - m ' T ^\ 



IA2I 2 ) 



det^/ 2 



S + M ~z 



det w / 2 



S + M I* 
~z* S^+M^ 



(21) 



where z — diag(z, — z) with z = (// + iT)l^,lN . The square root over the determinant is given because the number 
of the fermion measures is a half of that of the Nambu-Gorkov basis. We have defined the effective potential f2 as the 
function of the order parameters. In the thermodynamic limit, N — > oo, the ground state of the model is determined 
by the set of the order parameters which minimizes the effective potential. 

For the ground state solutions, we make a few Ansatzs for the order parameters. First, we set to be real, 



b*j = 4>f. Second, we also set A^ and A^ to be real and — A^ 



A« 



A a- Both assumptions are consistent for 



the formation of the scalar and parity-positive condensates in the ground state, which are favored by the finite quark 
mass term and the real r] A . Then, the effective potential is a function of six order parameters, <pf with / = u,d, and 
s and A A with A = 2, 5, and 7. Using these assumptions, the effective potential f2 becomes 

A 



D 



n =— 



3^ ' 3 
1 

8N c N f 



-rA 



[lndet(5 + M + z)+ lndet(S' t + M f - z) + Indet(5 + M + z*) + lndet(5 t + M ] - z*)] , (22) 



1/2 



7 



where A = 3E 2 /(2A^ c A^/Ga) and B = 3E 2 / \2N C N jG x ) are defined. We can obtain the ground state by solving the 
six gap equations, dtt/d(j)f = and dfl/dA^ = 0, simultaneously. 



D. Remarks 



To relate the order parameters in the ChRM model with the expectation values of the fermion bilinears in the 
microscopic theory, we use the external field derivatives as 

The order parameters 4>f and are proportional to the chiral and the diquark condensates respectively, and then 
we simply use the values of 0/ and A^ to distinct each phases. 

Note that, in the partition function, the parameter E can be absorbed by rescaling the order parameters, as well 
as the parameters, T, /j,,mf, and tja- Therefore, in the chiral limit rtif — (together with tja = 0), a change of E 
affects only on the scale of the phase diagram, and the global structure of the phase diagram is invariant. In fact, the 
parameter that can change the structure of the phase diagram is B/A, which is independent of E. In our treatment, 
the ratio is fixed by the Fierz coefficients, and is obtained as B/A = 3/8 with N c = 3. In Ref. various structures 
of the phase diagrams has been found with B/A changed. 



1 dh\Z _ 2B 

'2NN c N f dm f ~ (23) 

1 d\nZ 2A 
^ = — A4. 24 

2NN c Nr drj A 3 V ' 



IV. THREE EQUAL-MASS FLAVORS 

We first examine the ground state in the case with the exact flavor SU(3) symmetry. For this purpose, we set 
m u = rrid = 7Ti s = m. Assuming that the flavor symmetry is not broken spontaneously, we can set the order 
parameters as <f> u = 4>d = <f> s = 4> an d A2 = A5 = A7 = A. 

The effective potential is simplified to the function of the two order parameters, 

Cl = AA 2 + B(j) 2 - -t- ^2 ln[(<7 ± zf + A 2 ] 8 [(a ± zf + (2A) 2 ] + c.c, (25) 
72 ± 

where a = <j) + m and we have set r\A = 0. Combining two gap equations, dfl/d(f> = and dfl/dA = 0, we can 
determine a ground state solution for given T, ii 1 and m. 

It is easy to find that A = is always a solution of the gap equation, since A appears as A 2 in the effective 
potential. Moreover, if m = 0, <fi = is also a trivial solution for any T and /1. Then, in the chiral limit, there are 
generally four types of the solutions, (i) = 0, A = 0, (ii) <f> ^ 0, A = 0, (iii) (4 = 0,A/0, and (iv) ^0,A/0. 
When 777 7^ 0, the solution <p = no longer exists and is replaced by a small value proportional to m. 

Let us first consider solutions with A = 0. When A = 0, the effective potential (|2"5)) becomes identical with that 
analyzed in Ref. [i~T| . 

n = ^-\Ha 2 -z 2 )-\Ho 2 -z*% (26) 

where G 2 = 2B is defined. Therefore, the phase diagram described by the effective potential ([25]) at A = is the same 
as in Ref. (Tlj . Several points on the phase structure with the effective potential (|26|) are summarized: When m = 0, 
we find a second-order phase boundary in the large T and small \i region, and a first-order in the small T and large 

fi region. Two lines are connected at the TCP, (T 3 ,/i 3 ) = (^Vv^ + 1G~\ \\/ \/2 - 1G" 1 ) = (0.776G" 1 , 0.322G^ 1 ). 
The transition temperature at \i — is obtained as To = G _1 , while the transition chemical potential at T = is 
/Jo = 0.528G" 1 . We use these two values of To and /xo for a normalization of T and \x in the presentation of the phase 
diagram to remove the E dependence as possible. If finite m is introduced, the second-order phase transition line 
becomes smooth crossover, while the first-order line remains robustly. The TCP also becomes the critical point. 
We next consider the solutions with <j> = 0. By setting = and m = 0, the effective potential becomes 



£1 = AA 2 - - 
4 



^ln(A 2 + z 2 ) + iln(4A 2 + z 2 ) 



+ c.c. (27) 
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0.5 1 1.5 0.5 1 1.5 

n/n n/n 



FIG. 1. the phase diagrams with three equal-mass flavors. The left (right) panel shows the case with m — (Em = 0.1). 
The solid lines denote the first-order phase transitions, and the dotted lines the second-order phase transitions. The largest 
order parameters in each phases are shown by its letter. Tb (/j,q) is defined as the chiral phase transition temperature (chemical 
potential) on the T (fi) axis in the chiral limit when A = 0. 



The gap equation for a nontrivial solution A ^ is obtained as 

1 1 4 



9A 2 + z 2 9 4A 2 + z 2 



+ c.c. =0. (28) 



For a large T and/or /z, this equation does not have a real solution of A, which indicate that at some values of T and 
fj,, its solution coalesces to the trivial solution A = 0, where the system reaches a second-order phase transition. It is 
easy to find a curve of the phase transition by setting A = in Eq. 



" 2 - T2 - 2A (29) 



(fi 2 +T 2 ) 2 3 ■ 

The true phase structure should be, of course, determined by comparing the effective potential for all solutions of the 
gap equations, and the second-order phase transition line (|29p may be replaced by other phase structures. 

To investigate the whole phase diagram, we have to numerically compare the effective potentials for all possible 
solutions. The result is presented in Fig. [1] for the case with m = (left panel), and with ra/0 (right panel). 

In the chiral limit, we find a A = solution becomes the ground state for a small chemical potential [i/no ^ 0.95. In 
this region, the phase diagram is the same as that investigated in Ref. : There is the second-order phase transition 
line between the T axis and the TCP, and the first-order line from TCP to the large /j, region. 

When a chemical potential exceeds a critical value, however, we find a first-order phase transition to the CFL phase 
at a small temperature. In this phase, the chiral order parameter </> becomes zero. We stress that the phase transition 
to the CFL phase in the ChRM model is remarkable because in the construction of the ChRM model, the ratio of 
the couplings between the quark-antiquark and the quark-quark channels is given by the symmetry consistent with 
QCD, not tuned so that the phase transition can be reproduced. 

Unfortunately, on the other hand, the CFL condensate A continuously goes to zero not only as T is increased, but 
also as /i is. In QCD, the second-order phase transition at a large [i is not expected, since the Cooper instability 
remains even if the attractive interaction is infinitely small. This unphysical phase transition may be explained by the 
absence of the Fermi surface in the ChRM model, which is regarded as a model without the spacial dimension. We 
consider that, at such a region, the model reaches a limitation. Similar structures are found in the ChRM models for 
two-color QCD with the diquark baryon condensates [13, [H| > as wei l as in the study on the 2SC phase in the ChRM 
modelpj]. 

When m ^ (the right panel of Fig. [!}, we find qualitative and quantitative changes from the case with m = 0. 
Due to the nonzero m, <f> has a small value even in the symmetric and the CFL phases. The second-order chiral 
phase transition line is washed out to become a crossover, and then the TCP becomes a critical point. Note that the 
second-order phase transition line for A remains. The first-order chiral phase transition line is pulled up in the larger 
T and /i directions. The phase transition line between the chirally-broken phase and the CFL phase also shifts to 
the larger /i region. Although the second-order CFL phase transition line moves to extend the CFL phase, the CFL 
phase shrinks in total. 
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V. 2+1 FLAVORS 



We next concentrate on the response of the model to the asymmetry between the light up and down (ud) quark 
flavors and the mid-light strange quark flavor. In order to sec this effect, we set the quark masses as m u = rrid = m 
and m s 7^ m. Wc assume that the flavor symmetry is not broken further spontaneously, and set the order parameters 
as <t> u = §d — <t> an d A5 = A7 = A s . We also write A2 = A for convenience. Due to the asymmetry between the 
ud and strange quarks, the typical symmetry in the CFL phase SU c+ l + #(3) is not realized. Nevertheless, we call the 
phase with A^O and A s ^ the CFL phase. Also, the phase with A^O and A s = is defined as the 2SC phase 

Under the parametrization, the effective potential becomes the function of the four order parameters as 

SI =y(A 2 + 2A S ) + |(20 2 + 4> s ) 

y 2 E ln [{ ( ff ± z ) 2 + A2 } 3 i( ff ± ± *) + A '} 4 {(^ ± zf{{a ± z) 2 + A 2 ) + 4A 2 (( ( x ± z)(a s ±z) + A 2 )}] 
+ c.c, (30) 

where a = <f> + m, u s = 4> s + m s , and we set 7/^4 = 0. For the general case with <j) ^ 4> s and A 7^ A s , the determinant 
part under the logarithm becomes complicated. By setting m s — m, (j> s — <fi and A s = A, we recover the effective 
potential in the three equal- mass limit (|25|) . Another interesting limit is the 2SC phase, where A s = 0, in which the 
effective potential is separated to the ud quark sector and the strange quark sector as 

n = n ud (^,A) + n,(^) > (3 1 ) 

where 

Q ud (^, A) = | A 2 + ^0 2 - ^ ^ [4 In ((a ± z) 2 + A 2 ) + 2 ln (a ± z) 2 ] + c.c (32) 



n.OM = -fi - i Yl ln ± z f + c - c - ( 33 ) 



3 " s 72 



and 



The effective potential for the strange quark flavor (|33| is equivalent to the one of the conventional ChRM model 
without diquark condensates [TlJ, whose phase structure is summarized in sec. IV. 

The effective potential (|32|) can be compared to the one in Ref. [f4[ , where two light quark flavors are introduced and 
the strange quark degree of freedom is neglected. We first point out that the ratio of the coefficients of the quadratic 
terms of the order parameters, which can be read from (|32|) as (2B/3)/(A/3) = 2B/A = 3/4 since B/A = 3/8, is 
equal to that appearing in the model in Ref. [l4j |. The reproduction of 3/4 is important since the phase structure is 
sensitive to this ratio. 

Interestingly, however, these two effective potentials are not completely equivalent, and match only if T = or 
A = 0. This is because we use the different temperature scheme, or the matter effect matrix from that used in 
Ref. [3| ■ In the scheme used in Ref. [l4| , the effective temperature is introduced with opposite signs for two flavors^. 
On the other hand, in the scheme we used here, the temperature is introduced proportionally to the identity in the 
flavor space. The relation of these schemes is discussed in detail in Ref. [2l[. In Ref. (2f|, it is found that there is a 
unitary matrix in the flavor and the Matsubara space, which is a subspacc of the zero-mode space, to rotate the one 
scheme to the other. The fcrmion bilincars corresponding to the chiral condensates are invariant under the rotation, 
but the diquark condensates are not. This is the formal reason why we find that the two effective potentials are 
equivalent not only at T = 0, but also at A = 0. 

By considering the microscopic theory, a reason of the choice of the flavor-antisymmetric scheme can be explained 
as follows (2l|: To make the flavor-antisymmetric quark pair condensate be independent of the imaginary time r, only 
the terms of two quark fields with the opposite signs of the Matsubara frequencies should be nonzero in the Fourier 
summation over the Matsubara frequencies. This indicates that the temperature term, or the lowest Matsubara 
frequencies, should have the opposite signs for different flavors, because, in the diquark pairing, two quarks have 
different flavors in our treatment. 



4 Superficially, it seems to break the flavor symmetry, but, indeed, the symmetry is hold if the chemical potentials for two flavors are 
equal [2l. 
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FIG. 2. The phase diagram with two massless and one massive quark flavors, m u — = and Em a = 0.14. Notations are 
the same with Fig. 1. In the narrow slit between the symmetric phase and the CFL phase, which can be seen in the right panel, 
only A a has a large value. 



Nevertheless, we use the flavor-symmetric scheme in this paper because of the following reasons. First, if we 
construct the ChRM model relying only on the symmetry, we can not find a principle which scheme has to be 
adopted. Both treatments are consistent with the QCD symmetry, the anti-hcrmiticity at fj, = 0, and the chiral 
symmetry for all T and /j, of the Dirac operator. Second, because we treat the ChRM model in the case with three 
flavors, in order to adopt the flavor- antisymmetric scheme, we have to concern contributions from three combinations 
of the three flavors. This might make the effective potential more complicated, and we try to make the model as 
simple as possible. Indeed, the resulting phase diagram, which will be shown below, is qualitatively equivalent to 
the model with the flavor-antisymmetric scheme at A s = 0. In other words, they have the same global structure, or 
topology, of the phases. It suggests that the qualitative structure of the phase diagram is not sensitive to the selection 
of the temperature dependence schemes, as long as they hold the symmetry. 

We next present the resulting phase diagram. Similar to the case with three equal-mass flavors, we need numerical 
calculations to evaluate the whole phase diagram. The result with m = and m a ^ is shown in Fig. [5] At a small 
chemical potential, the diquark condensates do not have finite values in the ground state, and the phase diagram is 
the same with the model without the diquark condensates. Because for the ud quark flavors, there are no symmetry- 
breaking mass terms, the second-order phase transition line, as well as the TCP, exist. On the contrary, for the 
strange quark flavor, finite m s breaks the chiral symmetry explicitly, and then the second-order phase transition line 
becomes crossover. The TCP also becomes a critical point. In the symmetric phase, (j>, A and A s become exactly 
zero, but <j) s never becomes zero due to the finite symmetry breaking term m s . 

As n is increased with T kept small, a first-order phase transition to the 2SC phase is first observed. In the 2SC 
phase, not only A, but also <j> s have large values. If \x is further increased, we next observe a phase transition to the 
CFL phase, where both A and A s have large values. Note that because the flavor symmetry is explicitly broken, 
generally A ^ A s in the CFL phase. The diquark condensates become zero continuously when T and fi are increased. 
As seen in the right panel of Fig. [21 there is a narrow slit between two second-order phase transition lines for A and 
A s . In this area, only A s has a large value. Four phase transition lines, two horizontal second-order phase transition 
lines for A, and the vertical first-order phase transition lines for cj) s and A s , meet at one point. The second-order 
phase transition line for A s touches to the first-order phase transition line for <fi s . 

We consider the sequence of the melting of the diquark condensates. In the CFL phase, as T is increased with fi 
fixed, we find that the A becomes zero firstly, and the A s secondary. In the analysis of the NJL models 0, HH, on 
the contrary, it was found that A s melts firstly and A secondary, and then the phase with only A s finite is not found. 
The 2SC phase (where only A finite) above the CFL phase may also be expected in QCD. The reason is explained 
as follows]!, : Around the melting temperature, the quark Fermi spheres are smeared and then the sizes of gaps 
are mainly dominated by the density of states. Because the larger m s makes the strange quark Fermi sphere smaller, 
its density of states is smaller than that of the ud quarks, which results in the smaller A s than A. This indicates the 
preceding melting of A s . We consider that the contradicting result of our model is due to the absence of the Fermi 
sphere in our treatment, and, unfortunately, that this is a limitation of our model. 



11 



VI. SUMMARY AND DISCUSSION 



We have studied a ChRM model which can treat the competition between the chiral and the diquark condensates 
in the case with three flavors, to investigate the phase diagram with the CFL phase. In order to describe the color 
superconductivity, we have introduced the color and Lorentz indices to the truncated Dirac matrix. The random 
matrices mimicking the gauge fields are taken to be real matrices, whose elements are distributed according to the 
Gaussian weight. After the integration over the random matrices, we are left with the four fermion interaction term, 
which contains not only quark-antiquark interaction vertices, but also quark-quark vertices. Using Fierz transforma- 
tions, the ratio of the coefficients is uniquely determined. Applying the bosonization techniques, we finally derive the 
effective potential as the function of the order parameters, the chiral condensates and the diquark condensates. 

The phase diagram on the T-fx plane is calculated by solving the gap equations simultaneously. In the case with 
three equal-mass quark flavors, we find the CFL phase in the large chemical potential region, while the chirally-broken 
phase in the small chemical potential region. In the region where the CFL condensate is zero, the phase diagram 
is equivalent to that obtained from the conventional ChRM model without the diquark condensations. The CFL 
condensates become zero continuously, as T and/or /i are increased. The unphysical phase transition at large fi may 
be considered as the result of the fact that the ChRM model does not contain the Fermi surface. When finite quark 
mass m is introduced, the phase transition line for the CFL phase moves towards the larger chemical potential region. 
We also find that the region of the CFL phase becomes smaller as m is increased. 

For the case with 2+1 flavors, m u = md = and m s ^ 0, we find both the CFL and the 2SC phases. Moreover, we 
find the phase where A s has a finite value and A = 0. Such phase is not found in the NJL model, and not expected 
in QCD. We also consider this phase may be an model artifact due to the absence of the Fermi surface, but can not 
be excluded from the viewpoint of the symmetry. 

Although there are unphysical points, the ChRM model studied in this paper can consistently address both the 
chirally broken phase and the CFL phase. In addition, with the finite asymmetry between the ud quarks and the 
strange quark, the model also can show the 2SC phase. Because the ChRM model includes the Dirac matrix, 
investigation of its eigenvalue distribution, the Dirac spectrum, may be possible. It may be possible to understand the 
phase transition into the diquark-condensed states in the context of the moving of the Dirac eigenvalues. Furthermore, 
in the microscopic region, there mig ht be the universal structure of the Dirac spectrum, which can be compared to 
the other models at high density }24l 125}. The applications in these directions are postponed to future studies. 

Finally, we consider an outlook or a possible extension of the model. One of the most important effect we neglect 
is the U(l) breaking axial anomaly effect. It is known that, in the ChRM model, the anomaly effect introduces the 
flavor mixing, which changes the phase diagram drastically (iH l26l. |2~t|. Indeed, we find that the effective potential 
in the case with 2+1 flavors and A s = can be separated to the ud quark sector and the strange quark sector, 
as a result of the absence of the mixing effect. Moreover, on the phase diagram around mid chemical potential, it 
is suggested that the anomaly effect, which also mixes the chiral and the diquark order parameters, plays a crucial 
role on the phase structure with the various superconducting phases (28l-i30i]. We are then interested to combine the 
treatment of the ChRM model with the color superconducting phases presented here, and of the ChRM model with 
the axial anomaly in Ref. to investigate the effects of the axial anomaly on the phase structure with the color 
superconducting phases. It will allow us to discuss the phase structure at finite temperature and density from a 
viewpoint of the symmetry. 
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